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ABSTRACT 

We report on the detection of high-energy y-ray emission from the Moon during the first 24 months of observations 
by the Fermi Large Area Telescope (LAT). This emission comes from particle cascades produced by cosmic- 
ray (CR) nuclei and electrons interacting with the lunar surface. The differential spectrum of the Moon is soft 
and can be described as a log-parabolic function with an effective cutoff at 2-3 GeV, while the average integral 
flux measured with the LAT from the beginning of observations in 2008 August to the end of 2010 August is 
F{> 1 00 MeV) = (1.04 ± 0.01 [statistical error] ± 0.1 [systematic error]) x 10 -6 cm -2 s -1 . This flux is about a 
factor 2-3 higher than that observed between 1991 and 1994 by the EGRET experiment on board the Compton 
Gamma Ray Observatory, F(> 100 MeV) ~ 5 x 1 0 " 7 cm 2 s _1 , when solar activity was relatively high. The higher 
y-ray flux measured by Fermi is consistent with the deep solar minimum conditions during the first 24 months of 
the mission, which reduced effects of heliospheric modulation, and thus increased the heliospheric flux of Galactic 
CRs. A detailed comparison of the light curve with McMurdo Neutron Monitor rates suggests a correlation of the 
trends. The Moon and the Sun are so far the only known bright emitters of y- rays with fast celestial motion. Their 
paths across the sky are projected onto the Galactic center and high Galactic latitudes as well as onto other areas 
crowded with high-energy y-ray sources. Analysis of the lunar and solar emission may thus be important for studies 
of weak and transient sources near the ecliptic. 

Key words: astroparticle physics - gamma rays: general - Moon 
Online-only material: color figures 


1. INTRODUCTION 

The brightest object in the night sky, the Moon has been 
studied by naked eyes for millennia and is well known to shine 
by reflected sunlight. At hard X-ray and y-ray energies, the 
Moon shines as a result of cosmic-ray (CR) interactions with 
the lunar surface. Nuclear emission lines due to the spallation 
and activation of the surface material (regolith) by secondary 
neutrons as well as from natural radioactivity were first observed 
by y-ray spectrometers on Luna 10, 11, and 12 (Vinogradov 
et al. 1966, 1967) in the beginning of the space era and then by 
Apollo 15 and Apollo 16 missions (Metzger et al. 1973). These 
data provided the first insight into the chemical composition of 
the regolith. More recently, spatially resolved observations of the 
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nuclear emission lines by the Lunar Prospector (e.g., Lawrence 
et al. 1998; Prettyman et al. 2006), Kaguya (Kobayashi et al. 
2010), and Chang’E-1 (Zhu et al. 201 1) spacecraft were used to 
globally map the composition of the Moon. 

The high-energy continuum y-ray emission from the Moon 
was first detected by EGRET on the Compton Gamma Ray 
Observatory (which operated from 1 99 1 to 2000). Early analysis 
of the 1991-1994 EGRET observations yielded the integral flux 
of the lunar emission of F(E > 100 MeV) = (4.7 ± 0.7) x 
10” 7 cnT 2 s' 1 (Thompson et al. 1997), while a later reanalysis 
confirmed the detection and yielded a flux F(E > 100 MeV) = 
(5.55 ± 0.65) x 10” 7 cm” 2 s” 1 for the same viewing periods 
(Orlando & Strong 2008). 

The steady continuum y-ray emission of the solar system 
bodies arises from CR cascades developing in their surface 
layers. The principal emission processes are the decay of 
secondary neutral pions and kaons, produced by CRs hitting 
the surface, bremsstrahlung of primary and secondary electrons 
and positrons, and Compton scattering of the cascade photons. 
Although similar physical processes are involved for each body. 
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the observed y- ray spectra are not the same owing to the 
different environments for development of CR cascades, such as 
the gaseous atmospheres of the Earth (N, O) (Abdo et al. 2009) 
and the Sun (H, He) (Abdo et al. 2011), and the atmosphereless 
surfaces of the Moon and asteroids, y-ray emission lines have 
been observed from the rocky surfaces of Mars (Surkov et al. 
1992; Boynton et al. 2007), Eros (Evans et al. 2001), Vesta 
(Prettyman et al. 2011), and Mercury (Rhodes et al. 2011), 
though the Moon is the best studied so far. 

Calculations of interactions of CRs with the lunar surface are 
fairly straightforward, involving known physics processes, the 
well-measured spectrum and the composition of CRs near the 
Earth and composition of regolith, which was quite rigorously 
studied using the samples returned by the lunar missions as 
well as by remote sensing. The first calculation of the lunar 
y - ray emission including a continuum spectrum above a few 
MeV was done by Morris (1984) using the CR spectrum, 
cross section data, and techniques available at that time. Recent 
detailed calculations were made by Moskalenko & Porter (2007) 
using the framework of the Geant4 particle interaction code 
(Agostinelli et al. 2003). The simulations show that the spectrum 
of y-rays from the Moon is expected to be steep with an effective 
cutoff around 3-4 GeV (600 MeV for the inner part of the lunar 
disk) and exhibits a narrow pion-decay line at 67.5 MeV, perhaps 
unique in astrophysics. The emission from the limb of the Moon 
is expected to have a harder spectrum than that from the rest 
of the disk owing to the tangential directions of CR cascade 
development. Accurate measurements of the lunar emission can 
be used to test the details of the calculation with far better 
statistics and resolution than was possible with the previous 
EGRET observations. 

So far the Moon, the Sun, and the Earth are the only observed 
emitters of high-energy (> 100 MeV) y- rays in the solar system. 
Searches for other potential sources or populations of sources, 
such as small asteroids and comets, in the solar system are also 
ongoing. Having solid surfaces composed of rock or ice, such 
objects may exhibit a spectrum similar to the lunar spectrum 
at high energies, which can be a useful template. Populations 
of such small objects, which are too small to be resolved 
individually, could be detectable in the Main Asteroid Belt, 
Kuiper Belt, Oort Cloud, and near the orbits of Jupiter and 
Neptune (Jovian and Neptunian Trojans; Moskalenko et al. 
2008; Moskalenko & Porter 2009). 

As the Moon moves across the sky, the flux and spectrum of 
the Moon are needed for analysis of y- ray sources that it passes 
nearby. Although the motion of the Moon is fast compared to the 
integration times for most Large Area Telescope (LAT) source 
observations, it could affect analysis of flaring sources in close 
proximity. 

In this paper, we report observations of the Earth-facing 
side of the Moon by the Fermi - LAT during first 24 months of 
the mission. During this period solar activity was extremely 
low, resulting in a high heliospheric flux of Galactic CRs. 
Therefore, the CR- induced y- ray emission of the Moon was near 
its maximum. Preliminary results were reported by Giglietto 
(2009a, 2009b). 

2. OBSERVATIONS AND DATA SELECTION 

Fermi was launched on 2008 June 1 1 into circular Earth orbit 
with an altitude of 565 km and inclination of 25.6 and an orbital 
period of 96 minutes. The principal instrument on Fermi is the 
LAT (Atwood et al. 2009), a pair-production telescope with a 
large effective area (~6500 cm 2 on-axis at 1 GeV based on 


Table 1 

Summary of the Event Selection Cuts 


Cumulative Event Selections 

Photons 

(%) 

Live Time 

(%) 

Sroi ^ 20° 

100 

100 

\h i > 30 ° 

48.8 

56.6 

#Sun > 20° 

43.5 

50.2 

0(^1FGL > 2 x 10 -7 cm -2 s -1 ) > 5° 

40.1 

46.8 


post-launch Pass7 instrument response functions, IRFs) and 
field of view (2.4 sr or about 8 x 10 3 deg 2 ), sensitive to y- rays 
between 20 MeV and >300 GeV. After the commissioning 
phase, the LAT began routine science operations on 2008 August 
4. The Fermi - LAT normally operates in sky-survey mode where 
the whole sky is observed every two orbits. 

The analysis presented here uses post-launch P7V6 Source 
class IRFs. These take into account pileup and accidental 
coincidence effects in the detector subsystems that were not 
considered in the definition of the pre-launch IRFs. 711 The 
new IRFs also have improved effective area at low energies 
while the point-spread function (PSF) has been adjusted based 
on in-flight measurements. The systematic uncertainty of the 
effective area for the Source class events is estimated as 10% at 
100 MeV, decreasing to 5% at 560 MeV, and increasing to 20% 
at 10 GeV and above. The photon angular resolution is also 
energy dependent. The radius of the 68% containment region 
averaged over the LAT acceptance is about 4°2 at 100 MeV, 2:4 
at 200 MeV, 1°2 at 500 MeV, 0:67 at 1 GeV, and 0? 1 8 at 10 GeV. 

We use the LAT data collected between 2008 August 4 and 
2010 August 31, with a total duration of about 24 months cor- 
responding to solar minimum conditions. Our analysis of y- ray 
emission from the Moon is similar to the analysis of quies- 
cent solar emission described in Abdo et al. (2011). Events 
^100 MeV with a reconstructed direction within 0 20° from 

the position of the Moon (region of interest, ROI) and satisfy- 
ing the Source class selection (Atwood et al. 2009) are used. 
To reduce the contamination from the y- ray emission coming 
from CR interactions in Earth’s upper atmosphere, our selection 
is refined by selecting events with zenith angles <100° and ex- 
cludes time intervals when the ROI extends outside of this range. 
To reduce systematic uncertainties due to the bright diffuse 
y- ray emission from the Galactic plane and a possible spillover 
due to the broad PSF at low energies, we have also excluded 
the data taken when the Moon was within 30° of the plane (re- 
quiring |£>(j | ^ 30°). We further excluded the periods when 
the Moon was within 20° of the Sun (to exclude the brightest 
part of the extended emission due to the inverse Compton scat- 
tering of CR electrons off solar photons; Abdo et al. 2011) or 
within 5° of any other bright celestial source with integral flux 
Eifgl ^ 2 x I O' 7 cm 2 s -1 above 100 MeV as selected from 
the 1FGL Fermi - LAT source catalog (Abdo et al. 2010). These 
various selections produce a very clean event subsample, but at 
the expense of removing about 60% of the initial ROI data set 
as summarized in Table 1 . 

3. DATA ANALYSIS AND BACKGROUND EVALUATION 

Analysis of the Fermi - EAT observations of the Moon involves 
aspects not part of a standard point-source analysis. (1) The 
Moon moves across the sky rapidly, covering more than 10° per 


10 Details on the pre-launch instrument performance can be found in the 
Fermi - LAT calibration paper (Atwood et al. 2009). 
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day, and so the background for the source analysis is constantly 
changing. (2) The Moon is close enough to the Earth that 
parallax due to the motion of the Fermi spacecraft produces 
an apparent motion of the source relative to the fixed stars of up 
to 0. 6 in an orbit. (3) The Moon is not a point source. It subtends 
~0?5 on the sky. Such source extension, however, is measurable 
by the LAT only at energies above several GeV where the 
PSF is sufficiently small. At lower energies where most of the 
LAT photons are detected, a point-like source approximation is 
expected to be adequate. 

A dedicated set of tools was developed, not a part of 
the standard Fermi - LAT Science Tools package, to deal with 
moving sources such as the Moon and the Sun. Using these 
specialized tools, the data are selected in a moving frame 
centered on the instantaneous lunar position, which takes into 
account the parallax corrections and is computed using the JPL 
ephemeris libraries. 71 More details on the analysis of moving 
sources can be found in the paper by Abdo et al. (20 1 1 ) dedicated 
to the analysis of the y-ray emission from the quiet Sun. 

For the analysis of the Moon-centered maps, we used 
the Fermi - LAT Science Tools version 9r23pl. The gtlike 
tool provides maximum likelihood parameter values (using 
the method described in Cash 1979; Strong 1985; Mattox 
et al. 1996), which derives error estimates (and a full covari- 
ance matrix) from Minuit, a minimization tool supported by 
CERN, using the quadratic approximation around the maximum 
likelihood fit. 

As a first step, we verify the analysis procedure using simu- 
lated point-like sources with the computed spectra correspond- 
ing to (1) the inner part of the lunar disk (the total minus the 
limb emission), giving a softer spectrum, and to (2) the total 
flux, which has a harder spectrum (Moskalenko & Porter 2007). 
The analysis was done in Moon-centered coordinates using the 
same analysis chain and the same set of cuts as applied to the 
real data. Figure 1 shows the predicted differential y-ray flux 
for the inner part of the lunar disk and for the total emission 
corresponding to the solar minimum conditions together with 
the reconstructed spectra (points). The results demonstrate an 
excellent agreement between the reconstructed and simulated 
spectra. 

For a moving source evaluation of the background is of pri- 
mary importance. The main components of the background 
emission are the diffuse Galactic and isotropic (presumably ex- 
tragalactic) emission, weak point sources, and the instrumen- 
tal background (such as, e.g., misidentified CR particles). The 
evaluation of the background is made using the “fake-source” 
method described in detail in Abdo et al. (20 1 1 ). For this method, 
we apply to an imaginary source that follows the path of the 
Moon the same cuts as used for the Moon itself. Analysis of this 
“fake-Moon” yields an estimate of the background. To reduce 
statistical errors in the background determination, the back- 
ground is averaged over four fake-Moon sources displaced from 
each other and from the Moon itself by 30° intervals along the 
ecliptic. Because all fake-Moon sources are trailing the Moon 
at different angular distances over a period of 24 months, i.e., 
more than 24 full sidereal periods, they are sampling the same 
area on the sky. The background fluxes determined using indi- 
vidual fake-Moon sources are consistent within a fraction of a 
percent. Note that the 24 month analysis period is long enough 


71 http://iau-comm4.jpl.nasa.gov/access2ephs.html 
12 Available front Fermi Science Support Center (FSSC), 
http ://fermi .gsfc . nasa.gov/ssc 

73 http://seal.web.cern.ch/seal/work-packages/mathlibs/minuit/home.html 



Figure 1 . Verification of the analysis procedure for a moving source. Lines: 
predicted (Moskalenko & Porter 2007) differential y - ray flux for the inner part 
of the lunar disk (dashes, the total minus limb emission) and for the total emission 
(solid) corresponding to solar minimum conditions. Data points: reconstructed 
spectra from the simulated source. Open magenta circles — the central part of 
the disk and filled red points — total emission. 

(A color version of this figure is available in the online journal.) 

to average out any effects connected with incomplete sampling 
of the background. Attenuation of the diffuse background due 
to the shielding by the lunar disk itself is negligible compared 
to the observed flux because the solid angle subtended by the 
Moon is so small. 

4. RESULTS 

An ROI of 20° radius around the position of the Moon was 
selected, which is sufficiently wide to account for the broad 
LAT PSF at ~ 100 MeV (Atwood et al. 2009). Figure 2 shows the 
count profiles for events with E — 100 MeV-10 GeV within 10° 
of the Moon as a function of right ascension (a) and declination 
( S ) offsets with respect to the center of the lunar disk. It also 
shows the count distributions for the average fake-Moon and 
a sample of data simulated using a point-like source model. 
The simulation has been performed with the gtobssim simulator 
(LAT Science Tools), a standard tool for simulation of the 
y-ray sky. 

The angular size of the Moon is ~Ck 5. In order to test the 
possible effects of the finite angular size we plot the density 
of the events (100 MeV-10 GeV) around the Moon versus 
the angular distance from the center of the disk (Figure 3). 
Also shown are the event distributions found for the fake- 
Moon and simulated results for a point-like source with the 
spectrum matching the lunar spectrum. The event density has 
been computed as the number of observed events per annulus 
centered at the Moon position and normalized by the subtended 
solid angle. As expected, the angular size of the Moon has no 
effect on the angular profile of observed events due to the broad 
PSF of the lower-energy y- rays that dominate the data set. 

The distribution for a simulated source in Figures 2 and 3 
looks somewhat wider than the data, most likely due to devia- 
tions of the observed spectrum from the log-parabola fit (details 
of the fitting procedure are given below). Indeed, in Figure 4 
(left), the log-parabola fit seems to give a smaller number of 


74 Variations of the apparent size of the lunar disk due to its elliptical orbit 
average out over a time interval of about a month. 


4 


The Astrophysical Journal, 758:140 (8pp), 2012 October 20 


Abdo et al. 



projected onto the a (left) or S (right) offsets relative to the nominal position. Also shown are the count distribution observed for the fake-Moon and the simulation for 
a point-like source having the same observed spectrum. 



Figure 3. Event density vs. the angular distance from the center of the Moon 
for photons in the energy range 100MeV-10 GeV. The angular distance scale 
shows the inner radius of each annulus while the width of the annulus was fixed 
at 0?25 and the points are plotted at the mid-point radius of each annulus. Also 
shown are the count distributions observed for the fake-Moon and the simulation 
for a point-like source having the same observed spectrum. 


events than the observed spectrum at low energies. Since the 
integral flux of the simulated source matches the integral flux of 
the real source, this, combined with the energy dependence of 
the PSF, leads to shorter tails for the simulated distribution and 
somewhat higher flux at the core. The full maximum likelihood 
analysis of the data sample indicates that the maximum of the 
emission is centered at the computed position of the Moon, with 
a 68% confidence level error radius of 0(03. 

In order to derive the energy spectrum of the Moon, we deter- 
mined the background intensity using the previously described 
fake-Moon analysis, where we applied the same data selections 
that were used for the real Moon. The background was not 
allowed to vary in the likelihood analysis. Figure 4 shows the 
spectrum of the Moon, compared with past observations and the- 
oretical predictions. The observed spectrum was fit with a sim- 


ple power law and a log-parabolic function. The log-parabolic 
function 

dN/dE = N 0 (E/E b )-“-^ 1 °s (£/£ i. ) ! (i) 

provides the best test statistic 7 '' (TS). Here Nq = (70 ± 
6) x 10~ 9 cm' 2 s _1 MeV“\ a = 1.39 ± 0.07, ft = 0.45 ± 
0.02, and E b = 44.2 ± 2.3 MeV. This yields the integral 
flux F(>100 MeV) = (1.04 ± 0.01 [stat] ± 0.1 [syst]) x 
I O' 6 cm 2 s' 1 for a point-like source at the expected position, 
where the two uncertainties are statistical and systematic. Note 
that the parameters of the log-parabolic function and their errors 
are not independent; therefore, the relative statistical error of the 
integral flux is significantly smaller than the relative errors of 
individual parameters. 

The systematic error shown above has been evaluated by 
repeating the fit with the normalization of the background 
derived from the fake-Moon analysis left as a free parameter. 
It was found that in this case the normalization of the modeled 
background emission changes by <1%. The systematic error 
also includes the effect of uncertainties in the effective area at 
low energies. 

In addition to the spectral fit covering the full energy range, 
we constructed the spectrum for logarithmically spaced energy 
bins using the gtlike tool. For the fluxes in finite energy ranges, 
separate likelihood analyses were made for each range, with a 
power-law spectral model. Table 2 shows the differential flux, 
which was determined from the observed flux in the given energy 
range by dividing by the corresponding energy interval. The 
table also shows corresponding statistical errors and the TS 
values. Only an upper limit is shown for the entry with TS < 10. 

Since the y- ray emission of the Moon is induced by CR 
interactions, it should correlate with the ground neutron monitor 
(NM) rates, which measure the integral intensity of Galactic 
CRs. On the other hand, it should anti-correlate with the 
maximum extent of the heliospheric current sheet (tilt angle) 
that has been used in modulation models as a proxy for solar 
activity since the 1970s (e.g., Potgieter 2008). 

To test the correlation, the whole period of observations was 
divided into two-month intervals, which are long enough to 
keep the statistical error sufficiently low (Table 3). Figure 5 

75 The likelihood test statistic, TS = 2Alog(likelihood), between models with 
and without the source is defined in Mattox et al. (1996). 
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Figure 4. Left: y - ray spectrum of the Moon derived from the first 24 months of observations by the Fermi-LAT (red filled circles) and the fitted log-parabolic function 
Equation ( 1) (black dashed line). The blue squares show the EGRET data corresponding to periods of lower (open squares, EGRET Low) and higher (filled squares, 
EGRET High) solar activity, respectively (Thompson et al. 1997). The gray shaded region is the result of the reanalysis of the EGRET data (1991-1995) by Orlando 
& Strong (2008). Right: LAT spectrum compared with the predicted spectrum (Moskalenko & Porter 2007) corresponding to the period of low solar activity (solid 
line). The dotted and dashed blue lines show the separate spectra from the inner part of the lunar disk (0:48 diameter) and the limb, respectively, assuming an average 
diameter of 0(52. 

(A color version of this figure is available in the online journal.) 
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Figure 5. Integral y-ray flux (> 100 MeV) observed during the first 30 months (same selection cuts as in Figure 4) in two-month intervals compared with the McMurdo 
NM rates (left panel) and heliospheric current sheet tilt angles (right panel) for two Hoeksema models (e.g., Hoeksema 1992; Zhao & Hoeksema 1995), the classic L 
model (dashed line) and a newer R model (dotted line). L av and R av are the averages of the tilt angles calculated for the northern and southern hemispheres. Systematic 
errors in the photon flux are shown in black, and statistical errors are in red. The NM rates are corrected for the site air pressure. 

(A color version of this figure is available in the online journal.) 


Table 2 

The Observed Spectrum of the Moon 


Energy Interval 
(MeV) 

Differential Flux ± Stat ± Syst 
(cm -2 s -1 MeV -1 ) 

TS 

90-110 

(1.57 ±0.06 ±0.16) x 10 -8 

1577 

110-140 

(1.01 ±0.03 ±0.1) x 10 -8 

2746 

140-190 

(5.7 ± 0.1 ± 0.5) x 10 -9 

4576 

190-250 

(2.54 ± 0.07 ± 0.2) x 10 -9 

3469 

250-320 

(1.16 ±0.04 ±0.1) x 10 -9 

2356 

320-400 

(5.1 ±0.2 ±0.1) x 10 -1 ° 

1449 

400-550 

(2.2 ± 0.1 ± 0.1) x 10 -1 ° 

1395 

550-750 

(5.7 ± 0.4 ± 0.3) x 10 -11 

464 

750-1000 

(1.5 ± 0.2 ± 0.08) x 10 -11 

158 

1000-1250 

(5 ± 1 ± 0.3) x 10 -12 

49 

1250-1500 

(2.3 ± 0.7 ± 0.1) x 10 -12 

25 

1500-2000 

(8 ± 3 ± 0.5) x 10 -13 

14 

2000-3000 

u.l. 90% c.l. 3.4 x 10 -13 

<10 


(left) shows the observed lunar flux in each sub-period together 
with the two-month averaged rates from the McMurdo NM 76 
in Antarctica, which has a low geomagnetic cutoff. There is 
a tentative correlation of the trends, with the rates slowly 
rising from the beginning of observations until the end of 2009 
December. A ~5% drop in the NM rates in 2010 corresponds 
to a ~10%-15% drop in the observed lunar flux. Figure 5 
(right) shows the observed lunar flux in each sub-period together 
with the tilt angle 7 ' for two of Hoeksema’s models (e.g., 
Hoeksema 1992; Zhao & Hoeksema 1995). The upper and 
lower histograms represent his old (L model) and new model 
(R model), respectively. The difference between the models 
illustrates the systematic error associated with determination 
of the tilt angles. There is a hint of an anti-correlation of the 
observed lunar flux with the solar activity. These comparisons 
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Table 3 

The Observed Light Curve of the Moon 


Time Interval, MJD 

Integral Flux ( E > 100 MeV) ± Stat ± Syst 
(xlO- 7 cm -2 s -1 ) 

54682.7-54736.7 

9.9 ± 0.4 ± 1.0 

54740.0-54794.6 

11.1 ±0.4 ± 1.0 

54794.7-54850.7 

9.9 ± 0.4 ± 1.0 

54850.7-54901.2 

11.3 ± 0.4 ± 1.0 

54911.9-54960.4 

11.1 ± 0.5 ± 1.0 

54966.6-55014.8 

10.2 ± 0.4 ± 1.0 

55021.1-55069.4 

11.1 ±0.4 ± 1.0 

55075.6-55130.6 

10.7 ± 0.5 ± 1.0 

55130.7-55186.6 

11.0 ± 0.5 ± 1.0 

55186.7-55233.3 

10.2 ± 0.5 ± 1.0 

55243.0-55298.7 

10.7 ± 0.5 ± 1.0 

55298.7-55354.6 

10.2 ± 0.4 ± 1.0 

55354.7-55410.6 

9.6 ± 0.4 ± 1.0 

55410.7-55466.6 

9.1 ±0.5 ±0.9 

55472.3-55521.2 

8.5 ±0.4 ±0.9 

55521.2-55576.0 

9.7 ± 0.4 ± 1.0 


illustrate the expected dependence of the observed y - ray flux 
on the CR flux as measured by the NM, and on the solar activity 
as traced by the tilt angle, but the dynamic range in this study is 
too limited to make a definitive conclusion. 

5. DISCUSSION AND CONCLUSION 

For the first time the y - ray spectrum of the Moon is measured 
with high statistical precision between ~ 1 00 MeV and 2 GeV. 
The observed spectrum is somewhat softer than expected 
(Moskalenko & Porter 2007) and matches fairly well the 
predicted soft spectrum of the central part of the disk above 
200 MeV. Predicted spectra shown in Figure 4 (right) are 
from a new calculation that employs the code developed by 
Moskalenko & Porter (2007), but shows the emission from the 
inner part of the disk with the diameter of 0.48. This calculation 
was selected as the best match for the observed spectrum. The 
limb does not seem to give a significant contribution to the total 
flux, except below 200 MeV and above ~ 1 GeV, where the inner 
disk spectrum exhibits a cutoff. At low energies (<200 MeV), 
the calculated total flux agrees well with the observations. This 
is because the angular distribution of the CR-induced y-ray 
emission is broad at low energies and there is essentially no 
distinction between the inner disk and the limb (dashed and 
dotted lines). 

Although the LAT cannot separately measure the disk and 
limb spectra, we assume that the spectral discrepancy originates 
from an overestimate of the limb contribution. We therefore 
consider a possible deficiency in the modeling. In particular, the 
effect of the roughness of the lunar surface was not considered in 
earlier modeling, which assumed a smooth surface. The relevant 
size scale for y - ray production is the particle interaction length, 
which is ~90 g cm~ 2 for GeV protons in oxygen (the most 
abundant element in Moon rock). The latter corresponds to 
~50 cm for lunar regolith (density ~1.8 g cm " 1 ), implying 
that the distribution of intermediate-sized rocks and boulders is 
the most important factor. Attenuation by this surface debris of 
the flux of CR particles coming from tangential directions would 
result in a softer y-ray spectrum resembling the observations, 
while the emission produced by CRs with larger incident angles 
is less affected. 

Physical properties of the Moon surface have been a matter 
of practical interest since the glorious days of early space 


exploration and the beginning of robotic (Lunokhod 1 and 
Lunokhod 2) and human (the Apollo Lunar Roving Vehicles) 
transportation on the lunar surface (see, e.g., Heiken et al. 1991). 
Recent measurements by the Lunar Orbiter Laser Altimeter 
(Rosenburg et al. 2011) allowed a systematic evaluation of the 
lunar surface roughness to be done for the first time. It was 
found that most of the surface is characterized by a fractal-like 
roughness. This information can be used for a description of 
the surface characteristics for a future realistic modeling of CR 
interactions with the lunar surface. 

A comparison of Fermi-LKA and EGRET observations shows 
clear evidence of an anti-correlation of the lunar y-ray flux 
with the level of solar activity. Subdivision of the EGRET 
data into two samples, data taken during the solar maximum 
1991-1992 and data taken during the period of moderate activity 
1993-1994, indicated integral fluxes of 3.5 x I 0 -7 cm -2 s -1 and 
7 x 1(T 7 cm -2 s _1 , respectively, suggesting that the y- ray flux 
is higher when solar activity is lower (Thompson et al. 1997). 
The LAT observations came at a time of extremely low solar 
activity, and the 2008-2010 LAT flux is higher than any of the 
EGRET results, confirming the trend. 

During the period reported in this paper, the McMurdo NM 
rates were up to 50% higher than the rates reported during 
the EGRET mission (Thompson et al. 1997). Correspond- 
ingly, the analysis of the Fermi-LXT data yields for the in- 
tegral lunar flux F(> 100 MeV) = 1.04 ± 0.02 [stat] ± 
0.2 [syst] x 10 -6 cm -2 s _1 , about a factor of two increase 
over the EGRET flux averaged over all periods of observations 
(~5 x 1 CT 7 cm -2 s~ 1 ), in reasonable agreement with predictions 
for the CR environment associated with low solar activity. A 
comparison of the LAT light curve with the McMurdo NM rates 
shows a tentative correlation with the NM rates tracing the inte- 
gral flux of Galactic CRs and a provisional anti-correlation with 
the tilt angle of the heliospheric current sheet, which is a proxy 
for the solar activity. Future LAT observations of variations of 
the lunar y- ray flux over a significant part of the solar cycle are 
necessary to clearly establish the correlation/anti-correlation. 

The Fermi-LAT observations of the Moon have provided 
detailed flux and spectral information that can be used to 
compare with future, more detailed models. They can also serve 
as a useful template to search for emission from populations 
of small bodies in the solar system (Moskalenko et al. 2008) 
and are available for analysis of LAT sources that are close to 
the Moon during flaring activity. Although indirect, the lunar 
y- rays offer the advantage of monitoring the CR flux outside 
Earth’s geomagnetic field. 
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